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MAXIMUM  LIKELIHOOD  ESTIMATES  OF  LINEAR 
DYNAMIC  SYSTEM  PARAMETERS 


by 

T.  G.  RYALL 


SUMMARY 

The  parameters  of  a  linear  dynamic  system  are  estimated  by  using  the  maximum 
likelihood  method. 

.  Maximum  likelihood  estimates  are  asymptotically  unbiased  and  efficient,  that  Is  they 
are^good^ estimates.  Furthermore  the  maximum  likelihood  estimate  of  a  fmctien  of 
parameters  is  the  function  of  the  maximum  likelihood  estimates  of  those  parameters. 
Two  different  situations  are  studied: 

(/)  the  excitation  force  is  random; 

(2)  the  excitation  force  is  deterministic. 
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NOTATION 


Drawing  a  distinction  between  a  random  variable  and  a  possible  realisation  of  that  random 
variable  leads  to  a  cumbersome  notation. 

Contextual  information  should  make  it  clear  which  possibility  is  referred  to.  As  a  result,  in 
this  report  no  distinction  will  be  made. 

A  matrix  with  components  Mi} 

x  A  vector  with  components  x, 

E  The  mathematical  expectation  operator 

X  A  summation  over  all  observed  frequencies 

m 

n  A  product  over  all  observed  frequencies 

u*  The  complex  conjugate  of  ■ 

nT  The  transpose  of  n 

n  The  number  of  ensemble  averages 
m  The  number  of  discrete  frequencies  observed 


FORMAL  DEFINITIONS 


O^dm)  =  M  observation  of  output  at  frequency  cu. 

/u)(ht)  =  jth  observation  of  input  at  frequency  cu. 

Sample  estimate  of  average  output  power 

r.m  = 

Sample  estimate  of  average  input  power 

Sample  estimate  of  average  cross  power 

*.(*»)  =  ^  I(  0(V>/uV, 


Sample  estimate  of  coherence 


\ 


i, 


S 

■t  T 

•;  i 

::  1 


1.  INTRODUCTION 

This  report  is  concerned  with  the  maximum  likelihood  estimation  of  dynamic  system 
parameters.  The  particular  linear  systems  considered  are: 

(i)  an  elastic  aircraft  on  the  ground ; 

(ii)  an  elastic  aircraft  in  flight. 

(The  method  used,  however,  is  of  quite  general  application.)  The  aerodynamic  feedback  forces 
depend  in  a  linear  way  on  the  past  history  of  the  aircraft.  The  transfer  functions  for  the  “aircraft- 
air”  system  cannot  be  represented  by  rational  transfer  functions.  However,  rational  transfer 
functions  can  approximate  non-rational  transfer  functions  to  any  degree  of  accuracy  over  a 
limited  frequency  range.  Chapter  2  gives  a  brief  summary  of  the  properties  of  linear  systems 
pertinent  to  the  estimation  process.  Chapter  3  formulates  the  deterministic  likelihood  function 
for  excitation  for  single  and  multiple  channels. 

“Globally”  stable  numerical  algorithms  are  also  presented  to  solve  the  likelihood  equations 
and  an  estimate  is  made  of  the  asymptotic  covariance  matrix. 

Chapter  4  formulates  the  likelihood  function  for  random  excitation  of  single  and  multiple 
channels.  A  “globally”  stable  numerical  algorithm  is  presented  for  the  single  channel  case. 

No  attempt  is  made  to  solve  the  general  inverse  problem  of  determining  the  mass  density 
variation  throughout  the  structure  nor  the  boundary  conditions.  There  is  also  no  attempt  made 
at  finding  coefficients  of  a  set  of  linear  differential  equations  which  is  a  minimal  realisation  of  the 
process. 

The  purpose  of  the  estimation  procedure  is  to  provide  estimates  of  the  damped  modes  and 
frequencies  which  may  be  compared  with  the  results  of  finite  element  program  and/or  allow  the 
calculation  of  generalised  mass,  generalised  damping  and  generalised  stiffness.  The  Ho 
algorithm  (Anderson  and  Moore)  provides  a  method  to  obtain  a  minimal  realisation  from  a 
transfer  function  matrix.  It  is  well  known  that  the  general  inverse  problem  is  unstable,  i.e.  large 
changes  in  mass  density  variations  can  lead  to  small  changes  in  the  impulse  response.  An  extremely 
good  summary  of  the  maximum  likelihood  method  is  given  by  Kendall  and  Stuart  [I].  The 
asymptotic  distribution  of  the  finite  Fourier  transform,  on  which  this  method  of  estimation 
relies  for  robustness  is  covered  by  Brillinger  [2]  and  Anderson  [3].  Only  elementary  result  in 
dynamics  theory  are  called  for  with  the  exception  of  the  minimum-phase  transfer  functions  and 
the  fact  that  the  principal  subdeterminants  of  the  transfer  function  matrix  under  certain  con¬ 
ditions  form  a  minimum-phase  system. 

There  can  be  a  number  of  reasons  for  estimating  dynamic  system  parameters.  The  simple 
answer  being  to  obtain  a  mathematical  model,  however  this  really  begs  the  question.  Given  a 
mathematical  model  of  an  aircraft  on  the  ground  it  is  possible  to  compute  a  model  for  an  air¬ 
craft  in  flight.  Stability  questions  may  then  be  answered,  i.e.  for  what  range  of  values  of  forward 
speed  and  altitude  is  the  aircraft  stable?  Control  laws  may  also  be  devised  to  stabilise  an 
unstable  system  or  to  move  the  system  in  the  shortest  possible  time  (subject  to  control  constraints) 
to  another  state.  The  last  problem  requires  an  accurate  model  of  the  dynamical  system.  The 
literature  on  system  identification  is  immense;  a  good  survey  paper  is  System  Identification — A 
Survey,  by  K.  i.  Astrom  and  P.  Dykhoff.  The  paper  Spectrum  Analysts— A  Modem  Perspective  by 
S.  M.  Kay  and  S.  L.  Marple  Jr.  looks  at  some  of  the  deficiencies  of  various  estimation  procedures, 
especially  in  the  time  domain,  and  suggests  that  the  direction  of  future  research  is  in  formulating 
computationally  efficient  maximum  likelihood  auto-regression  moving  averages  for  spectral 
estimation.  This  report  is  confined  to  the  frequency  domain  and  as  such  has  a  number  of  inherent 
advantages  for  linear  systems: 

(i)  tightly  damped  modes  “separate”  more  than  in  the  time  domain; 

(ii)  measurement  noise  and  process  noise  are  much  easier  to  model; 


l 


a  1 
i 


L _ i 


i 

t 


(iii)  the  Fourier  tnntform  being  a  turn  of  a  large  number  of  sample  time  valuei  allows  the 
application  of  the  central  limit  theorem  to  give  a  Gaussian  distribution  for  both  the 
real  and  the  imaginary  components— in  other  words  there  is  a  statistical  robustness; 

(iv)  if  the  noise  processes  are  stationary  each  frequency  component  is  asymptotically  inde¬ 
pendently  distributed. 

The  disadvantages  of  frequency  domain  analysis  are: 

(i)  non-linearities  are  harder  to  deal  with  in  the  frequency  domain; 

(ii)  non-stationarity  of  noise  causes  frequency  components  to  be  statistically  correlated. 


2.  LINEAR  SYSTEMS 

In  this  section  a  brief  summary  will  be  given  of  the  essential  properties  of  linear  systems 
pertinent  to  the  estimation  process. 

Consider  a  second-order  linear  stable  system 

M*+gi+£Zs=®  (2.1) 

where  x  may  be  finite  dimensional  or  infinite  dimensional. 

If  an  impulsive  force  is  applied  at  the  fth  position  it  follows  that 
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fth  position 


(2.2) 


Equation  (2.2)  is  equivalent  to  (2.1)  for  t  >  0+  with  initial  conditions  x(0+)  =  0 
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fth  position. 


The  response  at  time  r  is 
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where  u‘s  and  s’s  satisfy 

[a’tf +«£+£] a  =  0. 

If  impulses  are  applied  at  all  other  positions  then  the  impulse  response  matrix  is, 
the  matrix  is  symmetric, 

£(0  =  £(m.V"+«: 

If  the  matrix  is  not  symmetric  then 

DM-  Id^V'+ettV') 

I- I 

where  ?  tetisfet 


l*2M+*&+&  r*  -  o. 


(2.3) 

(2.4) 
I  that 

(2.5) 

(2.6) 
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For  a  linear  system  with  general  input  1(() 

0(/)  =  t)I(t)<*  (2.7) 

where  0(1 )  is  the  reponse  vector.  On  taking  the  Fourier  transform  of  equation  (2.7)  it  follows 
that 


£<•/(!«)  (2-8) 

where  T(lu)  is  the  transfer  function  matrix.  T{U.t  is  also  known  as  the  Fourier  transform  of  the 
Green’s  function.  Equations  (2.5),  (2.6),  (2.7)  and  (2.8)  are  very  general  results  and  do  not  rely 
of  the  "lumped  parameter”  model  (2.1)  which  assumes  a  spatial  averaging  of  the  true  partial 
differential  equation  and  that  the  system  has  no  memory.  The  probability  density  functions  of 
observed  time  functions  become  much  simpler  in  the  frequency  domain  than  in  the  time  domain 
hence  all  calculations  are  done  in  the  frequency  domain.  The  parameters  to  be  estimated  are  the 
damped  modes  (modes  of  decay),  decay  rates  and  decay  frequencies.  A  number  of  situations 
are  studied: 

(i)  linear  system  excited  by  a  deterministic  signal ; 

(ii)  linear  system  excited  by  a  random  signal. 


3.  DETERMINISTIC  EXCITATION 
3.1  Detenaiaistk  Excitation,  Single  Inpnt  and  Single  Outpat 
In  this  Section  a  linear  dynamic  system  is  excited  by  a  deterministic  force  at  a  particular 
point  and  the  response  is  measured  at  the  same  or  another  point.  Typical  force  inputs  that  may  be 
used  are: 


(i)  F(l)  =  A  sin  ^,1 1  0  4KT; 

(ii)  F(t )  =  "Hammer  Wow”. 

Any  force  history  that  has  an  approximately  flat  power  spectrum  over  the  frequency  range  of 
interest  may  be  used.  The  relationship  between  measured  output  and  measured  input  is 

0,to,  =  (3-1-1) 

where  N(lm1  is  the  noise  at  frequency  a>  and  E(N{lm))  =  0. 

=  <r2(«>).  (3.1.2) 

The  noise  may  be  considered  to  be  made  up  of  three  components,  viz: 

(i)  output  measurement  noise; 

(ii)  effect  of  non-linearities  (non-linearities  are  treated  as  though  they  come  from  a  random 
noise  source), 

(iii)  process  noise. 

In  the  time  domain  nothing  can  be  said  about  the  probability  density  function  of  the  noise, 
without  making  some  very  strong  assumptions.  In  the  frequency  domain,  however,  it  is  assymp- 
totically  true  that  the  number  of  points  in  the  transform  increases: 

(i)  that  the  real  and  imaginary  parts  of  have  a  normal  probability  density  function; 

(ii)  each  frequency  component  is  independently  distributed. 

Statement  (i)  has  two  exceptions.  At  zero  frequency  and  at  the  NYQIST  frequency  the  imaginary 
part  of  iV(to)  is  zero  with  probability  one.  These  two  frequencies  will  be  excluded  from  this 
analysis  although  only  simple  modifications  are  needed  to  include  them.  The  likelihood  function 
may  be  formed  in  the  following  manner,  at  each  frequency  to  the  output  has  the  “complex” 
normal  probability  density  function 


I 


(exp-«0— r/XO*-f/*)/<r»}) 


where  <r2(w)  is  the  unknown  noise  power  at  to.  The  likelihood  function,  when  there  are  n  ensembles, 
is  at  a  particular  frequency  to 

(i)k4uj  (««P-mo-nxo*-r*n^c))) 

The  total  likelihood  function  for  n  ensembles  across  a  discrete  frequency  spectrum  is 

i-  =  n  (i)^J  exp-{l (O-T/KO*  -T*I*)lo\io)).  (3.1.3) 

The  likelihood  function  is  a  maiimmw  when  the  negative  log  likelihood  is  a  minimum.  Not  only 
does  L  contain  the  parameters  of  the  transfer  function,  L  also  contains  the  noise  spectrum  o\to). 

Now 


5ctj(oj) 

when 


a\to)  =  i  KO-TiyO-TI)*.  (3. 1 .4) 

n 

Substitution  of  (3.1.4)  into  (3.1.3)  and  taking  the  negative  logarithm  leads  to  the  minimisation  of 
•  (  " 

A  simple  interpretation  may  be  put  on  the  function  (3.1.S),  viz.  choose  parameters  in  the  transfer 
function  such  as  to  minimise  the  estimate  of  the  total  noise  in  decibels. 


(o-r/xo*-r*/») 


(3.1.5) 


3.2  Numerical  Solution  of  the  Objective  Faction 

In  any  optimisation  procedure  the  objective  function  should  be  well  scaled.  Let  the  para¬ 
meters  (decay  rates,  decay  frequencies  and  damped  modes)  of  interest  be  a,,  . . .,  up.  Then  the 
function  to  be  minimised  is 


f(«„  . . .,  *,)  =  i  £  logj^i  1(0-7/) .  (o-r/)*J 


where  m  is  the  number  of  frequency  points.  The  non-linear  equations  to  be  solved  are 


if 

6a, 


=  0 


1=  1 . P 


(3.2.1) 


The  Newton-Raphson  method  of  solving  non-linear  equations  states  that  if  a,  are  initial  estimates 
then  a,+Aa,  are  improved  estimates  where 


•f+X 

5a,  iajta. 


Aa,  =  0  for  all  f 


(3.2.2) 


Unfortunately  this  method  fails  to  converge  or  takes  a  long  time  to  converge  unless  the  initial 
estimates  are  sufficiently  close  to  the  true  values  so  as  to  ensure  that  the  Hessian  i*F  matrix 


is  positive  definite.  An  approximation  to  the  Hessian  matrix  is  used  which  is  globally  positive 
semi-definite.  This  approximation  is  the  asymptotic  expected  value  of  the  Hessian  matrix 


a F 
ia, 


hT  IT* 

_ zOl _ 

Mm  jjvp-muo-nr 


(3.2.3) 
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(3.2.4) 
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The  asymptotic  covariance  matrix  which  reflect*  the  amount  of  uncertainty  in  the  estimates  is 
given  by 


3.3  Detrrm  hlstic  Single  Pott  Exdtetien  MaWpk  Pelt  Mniwiatl 


If  a  linear  system  is  excited  at  a  single  point  and  measurements  are  mae  at  a  number  of 
points  then  there  are  two  possibilities  which  will  be  discussed: 

(i)  the  outputs  are  measured  simultaneously: 

(ii)  the  outputs  are  not  measured  simultaneously. 

If  the  outputs  are  measured  simultaneously  then  it  is  possible  that  the  system  output  is  con¬ 
taminated  by  correlated  noise.  Under  these  circumstances  the  maximum  likelihood  method  leads 
to  the  minimisation  of  the  “integral"  of  the  logarithm  of  the  determinant  of  the  estimated  covari¬ 
ance  matrix  of  the  noise.  The  function  to  be  minimised  in  the  case  of  two  outputs  it 


I  log 
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<0.-r../iX0»-rt,  /?)  (Oi-ruAXOS-rJ./T) 

(0a-r„/,x0?-r?./.j  (Oi-fj./.xoj-rj./T) 


(33.1) 


Sufficient  statistics  for  this  situation  are:  all  the  average  output  powers  and  all  the  average  cross- 
powers  between  outputs  and  all  the  average  cross-powers  between  output  and  input.  If  the  outputs 
are  not  measured  simultaneously  and  are  not  measured  during  the  same  excitation  then  the 
correlated  noise  is  assumed  not  to  be  present  and  the  off-diagonal  terms  in  (3.3.1)  vanish. 

The  maximum  likelihood  method  under  these  circumstances  leads  to  the  minimisation  of 
the  sum  of  all  the  “integrated"  logarithms  of  the  estimated  noise  power  spectra.  This  result  is 
the  direct  analogue  of  the  single  input-single  output  case. 


4.  RANDOM  EXCITATION 
4.1  Single  Output 

If  a  system  is  excited  by  random  noise  then  both  real  and  imaginary  components  of  the 
output  have  a  normal  distribution.  Each  frequency  component  is  asymptotically  independently 
distributed.  If  the  real  and  imaginary  components  of  the  input  are  independent  and  have  the  same 
expected  power  then  phase  has  a  circular  uniform  distribution  and  output  power  has  a  negative 
exponential  probability  density  function.  The  previous  restriction  is  equivalent  to  assuming  that 
the  process  has  a  certain  type  of  stationarity. 

Since 


0  =  TI 

(4.1.1) 

it  follows  that 

OO*  =  TT*/I* 

(4.1.2) 

Hence  the  likelihood  function  at  a  particular  frequency  is 

1 

u* 

(4.1.3) 

where 

7T*. 

(4.1.4) 
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If  there  are  N  ensembles  the  likelihood  function  is 


The  likelihood  function  across  a  discrete  frequency  spectrum  is 


/ l\*  f  zoo*  n 

\p)  e>  1  u*  py 

icrcte  frequency  spectrum  is 


Maximising  (4.1.6)  is  equivalent  to  minimising 
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(4.1.5) 

(4.1.6) 

(4.1.7) 


A  unique  solution  for  the  transfer  function  is  not  possible  even  if  P,a)  is  known  without  error. 
This  is  because  if  T(t)  =  N(JD(I)  where  N  and  D  are  polynomials  in  s  (the  Laplace  variable)  then 


PO) 


NQ)  N(—s) 
D(s)D(-s) 


(4.1.8) 


D(s )  is  uniquely  determined  since  the  system  is  stable.  NO)  is  not  uniquely  determined  since  any 
pair  of  conjugate  zeros  can  be  replaced  by  their  “unstable”  images  and  still  have  the  same  power 
spectrum.  If  there  are  4 N  complex  zeros  in  P(s)  this  means  that  there  are  2f2**)  possible  solutions 
for  TO)-  However  the  damping  and  frequencies  are  uniquely  determined  from  D(s).  As  an 
example  of  this  non-uniqueness  consider  the  two  time  histories : 


(i)  — e*3,+2e~4‘; 

(ii)  5«-J,-6e-*'. 

(4+ a>2) 

Both  these  time  histories  have  the  power  spectrum  density  Multiplication  of 

the  above  time  histories  by  —1  will  generate  another  two  time  histories  still  having  the  power 

44-tu2  . 

spectrum  density  - 5 — - 5- .  It  should  be  noticed  that  the  exponents  remain  the  same  in 

'  (9+<o  Xlh+«r) 

all  four  time  histories.  The  minimum  phase  system  is  -e"3'-t-2e~4'. 
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4.2  N— erica!  StMta  of  the  Likelihood  Equation 

After  appropriate  scaling  the  function  to  be  minimised  is 

|l[log  />,„>  +  ~  M  =  FOw,)  (42.1) 

06*  — — ■ 

where  M  is  the  number  of  frequency  points.  For  ease  in  notation  define  S(m)  =  where  OO* 

is  the  measured  averaged  output  power  and  II *  is  the  known  average  population  input  power. 
The  equations  to  be  satisfied  are 


4F 

da, 


i  —  I, .... p 


If  #|  are  initial  estimates  then  «,+ Aa,  are  improved  estimates  where 

dF  „  d*F  A 

—  +  £  - Aa,  =  0 

da,  daja ,  1 


(4.2.3) 


Now 


dF  _  Mt  ¥  _ 

>a,  ”  ujpda,  p2da,f' 
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(4.2.4) 
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i2F  = 

Sa^ay  '  A7(/>2  3a,  4a,] 


The  Hessian  matrix  is  again  approximated  by  a  globally  positive  definite  matrix.  The  matrix  is 
replaced  by  its  expected  value,  and  the  approximation  is  asymptotically  exact. 


System  Excited  by  aa  Unknown 


Prepare  FieM 


If  a  linear  system  is  excited  by  an  unknown  random  pressure  field,  not  necessarily  white  in 
space  and  time,  then  it  is  possible  to  recover  the  damped  modes,  frequencies  and  dampings  of  the 
system  provided  that  the  pressure  distribution  is  assumed  to  be  due  to  the  one  random  source. 
In  other  words  the  pressure  distribution  is  coherent.  In  theory  it  is  only  necessary  to  assume  that 
any  extraneous  noise  is  uncorrelated  between  measuring  stations.  However  to  avoid  working 
with  an  unmanageable  number  of  parameters  it  is  assumed  that  the  extraneous  noise  has  a  popu¬ 
lation  white  power  spectrum  and  this  is  of  the  same  magnitude  for  all  measuring  stations.  Let 

Oj  —  TJ-VNi  i  —  1, . . .,  m 

where  O,  is  the  measured  output  at  position  i  due  to  a  random  pressure  distribution  7  and  noise 
N,.  T,  is  an  overall  transfer  function  for  the  output  at  position  i  and 

E(N,N*)  =  a (4.3.1) 

£(//•)  =  I  E{lN*t)  =  0  (4.3.2) 

If  £(//*)  #  1  then  this  is  absorbed  in  T,.  Hence 

E(0,0*)  =  T,Tj+c%t.  (4.3.3) 

To  uncorrelate  the  data  it  is  necessary  to  find  the  eigenvectors  of  the  matrix  (4.3.3).  The  matrix 
(4.3.3)  has  a  degeneracy  of  order  m—  1  and  as  a  result  the  eigenvectors  may  be  defined  in  a  number 
of  ways.  This  is  now  done. 

Define 

Z,  =  IT?0,  (4.3.4) 

Zm  =  T,Om-TmO,  M>2  (4.3.5) 

then 

SCZ.ZJ)  =  (Zr.nXZT’.rT  +  e2)  (4.3.6) 

BKZ^Z)  =  (TiT\+T,T*.)o1  (4.3.7) 

EiZtf)  =  0  i+j. 

The  Z  random  variables  are  uncorrelated.  The  maximum  likelihood  method  in  this  situation  leads 
to  the  minimisation  of 

xjln(Er)rf)+ln(Xr,7'?-|-<72)-l-(/»— l)lna2+ 

£j[h*r1rT+r.r:)  +  +  (ir(r?xsr(rr+<r2)}- 

It  should  be  noted  that  the  ambiguity  in  the  transfer  functions  is  resolved  by  using  the  cross- 
power  spectra  (see  Appendix  A). 


S.  CONCLUSION 

The  maximum  likelihood  method  has  been  applied  to  the  estimation  of  linear  dynamic 
system  parameters  for  both  deterministic  and  random  excitation. 

Globally  stable  numerical  methods  have  been  presented  for  the  estimation  of  these  para¬ 
meters  together  with  an  estimate  of  the  corresponding  asymptotic  variance.  These  estimatse  have 
the  desirable  property  of  being  both  asymptotically  unbiased  and  asymptotically  efficient. 


APPENDIX  A 


Ambiguity  in  Transfer  Function  Resolved  Using  Cross-power  Spectra 

Let  T,(s)  and  T2(s)  be  two  transfer  functions  from  a  multiple  output  and  multiple  input 
model.  Let  /  be  the  random  input  to  each  (this  is  assumed  to  be  white).  Hence, 

O,  =  r,/  02  =  TJ  (Al) 

£■(0,0?)  =  r,rf£(//»)  =  r.rj  (A2) 

E(020\)  =  T.TtEan  =  r2rj.  (A3) 

T ,  and  T2  will  have  common  poles  and  these  may  be  identified  from  either  the  output  power 
spectrum  of  channel  1  or  the  output  power  spectrum  of  channel  2  using  the  knowledge  that  the 
system  is  stable.  Given  the  output  power  spectrum  the  zeros  of  TT*  may  be  determined  but  there 
is  a  question  as  to  how  to  assign  these  zeros  to  T  and  T*  respectively. 


r,  =  NtID  and  T2  =  N2ID 

where  Nt,  N2  and  D  are  polynomials  in  5.  It  follows  that 

E(020\)  =  n2nVd2. 


Hence 


(A4) 


E(020*)IE(020*2 )  =  Nj/Nf.  (A5) 

It  follows  that  provided  7\  and  T2  have  no  zeros  in  common  then  Tt  and  T2  may  be  determined 
from  power  spectra  and  cross-power  spectra.  In  practice  the  division  (A5)  is  not  performed  as  the 
argument  may  proceed  from  common  factors  of  E{020*)  and  E(02O\).  Overall  ambiguity  in 
sign  remains  since  multiplying  each  transfer  function  by  —1  does  not  change  the  power  spectra. 


L 


APPENDIX  B 


Sctatet  M—lllw  CoriiM  hi  a  Transfer  Faction  Matrix 

Consider  a  linear  time  invariant  vibrational  system  obeying  Maxwell's  law  of  reciprocity 

with  no  external  energy  source.  Let  T,j(s)  (i  =  1 ,  N;J  =  1 . N)  be  a  transfer  function  matrix 

then  the  following  statements  are  true. 

(i)  Each  element  has  the  same  poles  which  lie  on  the  LHS  of  the  complex  plane. 

(ii)  All  the  zeros  of  Tu  lie  on  the  LHS  of  the  complex  plane,  that  is  they  are  minimum  phase 
transfer  functions. 

(iii)  Any  principal  sub-determinant  forms  a  minimum  phase  transfer  function. 

(iv)  Knowledge  of  the  zeros  along  the  main  diagonal  determines  all  the  zeros  of  all  the  off- 
diagonal  elements. 

Comments  will  now  be  made  on  each  of  these  statements. 

(i)  The  location  of  poles  gives  information  on  the  damping  and  frequency  of  modes. 

It  is  possible  for  zero-pole  cancellation,  if  this  occurs  then  for  our  purposes  the 
zero  and  pole  are  both  replaced. 

(ii)  The  zeros  of  T„  correspond  to  the  dampings  and  frequencies  of  the  original  elastic 
system  rigidly  constrained  at  position  i.  Furthermore  these  frequencies  must  be  greater 
than  or  equal  to  the  frequencies  of  the  unconstrained  system,  for  lightly  damped 
systems. 

(iii)  The  zeros  of  any  principal  sub-determinant  correspond  to  the  original  elastic  system 
constrained  at  every  position  whose  element  is  left  along  the  main  diagonal.  As  the 
determinant  increases  in  size  the  transfer  function  tends  to  an  all-pole  filter. 

(iv)  As  an  example  consider  the  2  x  2  determinant 

T  T 

I  "  "  let  T,J  =  N,jlD  (r12  =  r21) 

\Tt»  T21 

where  N  and  D  are  polynomials  in  s.  The  determinant  is 
(NtlNi2-N\2)ID\s)  =  |£j 

hence 

N\2(s)  =  NtiN22-E(s)D(s) 

and  N\2  may  be  uniquely  determined,  provided  A,,,  N22  and  D  are  known.  £(s)  is 
determined  from  the  fact  that  NuN22~ED  has  repeated  roots  and 

N\i(S,)  =  N^SdN^S,) 

where  S,  is  any  zero  of  D{s).  The  zeros  of  any  off-diagonal  term  (TtJ)  are  calculated  by 
considering  the  binary  system 

j  Tij  | 

i  I 

\T,j  Tjj  \ 

The  total  transfer  function  matrix  cannot,  however,  be  reconstructed  due  to  ambiguity 
in  signs.  In  the  example  given  N2tJ  may  be  uniquely  found,  Ni2,  however,  has  an 
ambiguity  of  sign. 


f 


The  number  of  possibilities  considering  only  binary  systems  is  2(5).  Taking  into  account 
higher  order  determinants  reduces  this  number  to  2*“  *.  As  an  example  of  these  “phase’*  possi¬ 
bilities  consider  a  matrix  which  is  3  x  3  with  an  acceptable  phase  distribution 

'  +  +  +" 

+  +  + 

,+  +  +. 

then  the  other  possibilities  are 


'+  - 

— 1 

+ 

1 

4- 

_ i 

i 

+ 

+ 

1 

-  +  + 

-  +  - 

+  +  - 

_  —  +  +_ 

i 

+ 

1 

+ 

_ i 

“  +. 

This  pattern  of  signs  leaves  all  principal  sub-determinants  invariant. 
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